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Abstract 



Existence of homoclinic orbits in the cubic nonlinear Schrodinger equation under 
singular perturbations is proved. Emphasis is placed upon the regularity of the semi- 
group e^**^^ at e = 0. This article is a substantial generalization of and motivated 
by the effort of Dr. Zeng |^]. The mistake of Zeng in is corrected with a normal 
form transform approach. Both one and two unstable modes cases are investigated. 
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1 Introduction 



Consider the singularly perturbed nonlinear Sclirodinger equation, 

iqt = Qxx + 2[|gp - uj'^]q + ielq^,:, - aq + p] , (1.1) 

where q = q{t, x) is a complex-valued function of the two real variables t and x, t represents 
time, and x represents space. q(t, x) is subject to periodic boundary condition of period 
27T, and even constraint, i.e. 

q{t, X + 2tc) = q{t, x) , g(t, -x) = q{t, x) . 

C(j is a positive constant, a > and f3 > are constants, and e > is the perturbation 
parameter. 

In this work, we revisit the problem on the existence of homoclinic orbits in perturbed 
nonlinear Schrodinger equations 0. The crucial new feature is as follows: Singular pertur- 
bation ed^q will be investigated in contrast to the regular perturbation ed^q considered in 
0], where is a pseudo-differential operator obtained from a truncation of d"^. This study 
is motivated by the effort of Dr. Zeng 0] ^j. The main difficulty introduced by the singular 
perturbation ed"^ is that it breaks the spectral gap condition of the unperturbed system. 
Therefore, standard invariant manifold results will not apply. Nevertheless, it turns out 
that certain invariant manifold results do hold. The regularity of such invariant manifolds 
at e = is controled by the regularity of e^'^^ at e = 0. Difficulties and interesting results 
created by the singular perturbation term ed^q will all be commented in Remarks. 

The entire theory of locating a homoclinic orbit is divided into two parts. Part 1 
deals with local invariant manifold theory. Part 2 deals with global theory which includes 
integrable theory, Melnikov analysis, etc.. 

The notation | | will denote absolute value, and the notation || \\s will denote the Sobolev 

(i.e. W^'"^) norm of periodic function with period 27r. 

This article is written for experts only. Standard details will be omitted. 
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2 Local Theory 



Local theory is referred to a theory in a neighborhood of certain circle of fixed points, 
which includes local unstable fiber theorem, local center-stable manifold theorem, and size 
estimate of local stable manifold for certain saddle. These are some of the tools needed in 
locating a homoclinic orbit. 



2.1 Dynamics in a 2D Invariant Subspace 

The 2D subspace 11, 

n = {g| 9,^ = 0}, (2.1) 
is an invariant subspace under the PNLS flow { \L . 1|) . The governing equation in 11 is 

iq = 2[\q\^ - uj^]q + ie[-aq + p], (2.2) 



where ■ = 4i ■ Dynamics of this equation is shown in Figure |2ri|. Interesting dynamics is 



created through resonance in the neighborhood of the circle S'^^: 

S^ = {qeli \ \q\ = uj}. (2.3) 

When e = 0, S'oj consists of fixed points. To explore the dynamics in this neighborhood 
better, one can make a series of changes of coordinates. Let q = a/Jc*^, then (|2.2| ) can be 
rewritten as 

/ = e(-2a/ + 2pVlcose) , (2.4) 
e = -2(1 - u') - e(3'^ . (2.5) 

There are three fixed points: 

1. The focus Oe in the neighborhood of the origin, 



cos9 = ^, ^e(o,f). 



(2.6) 



P 

Its eigenvalues are 



sin 6* — ea, (2.7) 

where / and 9 are given in (P]B|). 



2. The focus in the neighborhood of S^i (|2.3|), 

cos^ = ^^, 0G(-f,O). 



(2i 
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Figure 2.1: Dynamics on the invariant plane 11. 



Its eigenvalues are 



/xi^2 = ±^v^\/-4y7/3sing + e ( I - ea, 



where I and 6 are given in ( |2.8|) . 
3. The saddle Qe in the neighborhood of S^^ (p.3|), 



cos^ = ^ 
Its eigenvalues are 



cos9 = ^, ^e(0,i) 



/xi,2 = ±v^W 477/3 sin^-e(^y 



— ea, 



where / and ^ are given in ( p.lO|) . 
Now focus our attention to order ^/e neighborhood of S^. 

J = I — uj"^, J = T = \ftt, 

we have 



uj^ + \fti cos d 



9' = -2j - V~ef3- 



sin 6 



where ' = ^ . To leading order, we get 
0' = -2j . 



3D and let 



(2.9) 



(2.10) 



(2.11) 



(2.12) 
(2.13) 



(2.14) 
(2.15) 



There are two fixed points which are the counterparts of Pe and (|2T8| ) and (|2.10 



Figure 2.2: The fish-like dynamics in the neighborhood of the resonant circle S^, 



1. The center P^, 

j = 0, cos6' 



TT 



ee (-2>o)- 



Its eigenvalues are 

/ii,2 = ±i2v^(/?2 - a^uj'^fi. 



2. The saddle Q*, 

i = 0, cos 9 



9e 



TT 



Its eigenvalues are 

fii,2 = ±2v^(/32-aV)3. 

In fact, (|2.14| ) and (|2.15|) form a Hamiltonian system with the Hamiltonian 

n = f + 2u{-auj9 + (3 sin 9). 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



Connecting to Q^, is a fish-like singular level set of 7i, which intersects the axis j = at 
and Q = (0,9), 



au{9-9,) = (3{sm9 -sm9,), 9 G (-y.O), 



(2.21) 



where 9^, is given in ( p.l8| ). See Figure for an illustration of the dynamics of ( p.l2| )- 
( ^.15| ). For later use, we define a piece of each of the stable and unstable manifolds of 

J = €{0), J = €{0), 9e[9 + 6,9. + 27T], 

for some small 6 > 0, and 
9 — 9^, 



<pm = -c(^^). 



■^/2uj[auj{9 -9^) -/?(sin^-sin^,)]. 



(2.22) 



0^(6*) and (pl{9) perturb smoothly in 9 and y/e into 0^ and for (|2.12|) and ( p.l3| ). 

The homoclinic orbit to be located will take off from along its unstable curve, flies 
away from and returns to 11, lands near the stable curve of Qe and approaches spirally. 
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2.2 Change of Coordinates 



As mentioned above, interesting dynamics happens in the neighborhood of the circle S^^ 
It is natural and convenient to center our coordinates around S^^- First, write q as 



g(t,x) = [p(t) + /(t,x)K^W, 



(2.23) 



where p and 6 are polar coordinates on 11 ( |2.1| ), and / has zero spatial mean. We use the 
notation (■) to denote spatial mean. 



27T 



2lT 



qdx. 



(2.24) 



Since the L^-norm is an action variable when e = 0, it is more convenient to replace p by: 



I={\q?)=p' + {\f?)- 
Since corresponds to / = cu^, the final pick is 



In terms of the new variables (J, 6, /), Equation (|1 . 1|) can be rewritten as 



J 



2a( J + cu^) + 2/3VJ + cj2 cos 6 + eUi, 
sin 9 



-2 J - e/3 



ft = LJ + VJ - iU2 - lAfs, 



where 



sin 6 



LJ- 

VJ = -i2JU + /) + ie[3f 



7^ 



-2(|/,|2) + 2/3 cos^ [v/J + a;2-(|/|2) _ ^/JT^ 
P 



eP sin 9 



Af, = 2p[2i\f\' - {im) + if - if))], 

^f3 = -{f+ p + 6i/r)/ + 2(i/iv - (i/r/)) 

--(l/P(/ + /))/-2(|/P)/ 



e/3 sin 9 



1 



(2.25) 
(2.26) 

(2.27) 

(2.28) 
(2.29) 

(2.30) 
(2.31) 

(2.32) 

(2.33) 
(2.34) 

(2.35) 
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Remark 2.1. The singular perturbation term "ed^q" can be seen at two locations, and 



TZ'l ( \2.3(\ , ^73^ ). The singular perturbation term {\fx\ ) in TZ^ does not create any difficulty. 



Since is a Banach algebra this term is still of quadratic order, (l/xP) ~ C^dl/lli)- 
Lemma 2.1. The nonlinear terms have the orders: 

\ni\-om% \K\-om% 

\W2\\s-omi). \\ms-om% (^>i)- 

Proof. The proof is an easy direct verification. □ 
2.3 Normal Form Transformation 

In locating a homoclinic orbit to (|2.10|) , we need to estimate the size of the local stable 
manifold of Q^. The size of the variable J is of order 0{^). The size of the variable 9 is 
of order 0{1). To be able to track a homoclinic orbit, we need the size of the variable / 
to be of order 0{e^), fi < 1. Such an estimate can be achieved, if the quadratic term N2 
( p.34| ) in (|2.29|) can be removed through a normal form transformation. 



In terms of Fourier transforms, 
and the two terms in A/2 can be written as, 

f - if) = E 



k+e^o 

1 1^ . . (2.36) 

^ E imf{-i) + f{i)f{-k)y^'^'>. 



2 

k+e^o 



It turns out to be convenient to work with the symmetrized form (|2.36| ). We will search 
for a normal form transformation of the general form, 

g = f + K{f,f), (2.37) 

where 



K{fJ)= E [Ki{k,£)f{k)f{£) + K,{k,e)f{k)f{-e) 
k+e^o 



+K,ii,k)fi-k)fii) + K,ik,i)fi-k)fi-i) 



e 



Kj{k,£), {j = 1,2,3) are the unknown coefficients to be determined, and Kj{k, 
K,ii,k), (j = l,3). 



9 



Lemma 2.2. For oj G (|, |) /S, S is a finite subset, there exists a normal form transfor- 
mation of the form ( p.37|) that transforms the equation 

ft = LJ-iAf2, 
into an equation with a cubic nonlinearity 

gt = L,g + 0{\\g\\l), > 1), 
where is given in ( p^.30D , and N'2 has the expression (cf: ( |2.34| ) ), 

Ar, = 2^[2(|/p-(|/r)) + (f -(f))]. (2.38) 
Proof. Denote the operator idt — iL^ by C^. We have 

C,g = CJ + CMf,f) 

= CJ - iL,K{f, f) + iKiLJ, f) + tK{f, LJ) 

+ iK{dtf - L J, /) + iK{f, dj - LJ), 

where Cef = A/2 and K{dtf — L^f, f) and K{f, dtf — L^f) will be shown to be cubic in /. 
To eliminate the quadratic terms, we need to set 

zL,K{f, f) - iK{LJ, f) - zK{f, LJ) = U2, 

which takes the explicit form: 



{ai+i(r)Ki{kJ) + bK2{k,i) + bK2{i,k) + bK3{-k,-i) = -2u, (2.39) 
-bKi{k, e) + {cx2 + ia)K2{k, £) + fei^2(-^, -k) + bKsik, £) = -2uj, (2.40) 
-bki{k, i) + bK2{-k, -i) + ((T3 + ia)k2{i, k) + bk^ik, i) = -2iJ, (2.41) 



where 



bKi{-k, -i) - bK2ik, i) - bK2ii, k) + {a^ + ia)K^{k, i) = 0, (2.42) 

6 = -2cJ^ a = e{2k£-a), ai = 2{ki + uj^), a2 = 2{f + ki - u^), 
as = 2{k^ + k£- 00^), a4 = 2{k^ + + k£ - 3uj^). 

Since these coefficients are even in {k,i), we will search for even solutions, i.e. 

k,i-k,-i) = k,ik,i), J = 1,2,3. 

(g)+(|^), (|^+(|^-, and (|^) + 0)- lead to 



{ai+ia)Ki{kJ) + bKi{k,i) = -K, 



((T2 + tcj)K2ik, i) - bK2ik, i) = -K, 



(as + ia)K2ii, k) - bK2{t, k) = -K, 
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where 



K = 2cu + {ai + i(j)k^{k, £) + hk^{k, i). 

Therefore we can express Kj{k,i) in terms of K as, 

ki{k,i) = {al + - b^)-^[bK - {cxi - ta)K], 

k2{k,i) = {al + or^ - \?)-\-hK - (as - %(y)K\, 

k^ii, k) = (al + a^- b^)-'[-bK - {as - ta)Ki 

ksikj) = {al + - b^)-^[{a^ - ia){K - 2uj) - b{K - 2u)]. 

Substituting these expressions into (|2.42|) , we get 

K = {\U\^ - \V\Y\WU -WV), 

where 



U 



V 



b{a2 — ia) b{a3 — ia) aj + a"^ 



+ 



+ 



af + a'^-b^ al + a^-b^ al + a^-b^ aj + a^ - b^ 



b{ai + ia) 



a; 



a^ 



52 



+ 



W = 2uj[a\ + a'^ 



^2 



a^ 



52 



+ 



an 



a^ 



62 



6(o"4 + ia) 
al + a^-b^' 



[0-4 + 0- - 6(0-4 + ia) 



(2.43) 
(2.44) 
(2.45) 
(2.46) 

(2.47) 

(2.48) 

(2.49) 
(2.50) 



For uj e (|,|), the denominators in (|2^ ) - (|2^) and (|2^ ) - (l2^) , and aj (1 < j < 4) 
vanish at u; in a finite subset. For u not in this finite subset, a"^ is a O(e^) small perturbation 
of o-| — 6^ (1 < j < 4), and a is a (9(e) small perturbation of aj (1 < j < 4). Setting a = 0, 
we have K = 2uj, which leads to the solution given in [R]. To the order 0{e), 



K = 2u 



1 + iba 
b 



1 



+ 



al-h' 
b 



+ 



al 



62 



+ 



+ 



(74 



0-2 



0-3 + 6 0-4 — 6 



62 
-1" 



(2.51) 



We will show that the denominator in ( |2.51|) does not vanish except for in a finite subset. 
Denote the denominator byZ). As/c— >±oo, or£-^±cxD, 



D ^l. 
We also know that 

D = l + b 



(2.52) 



1111 
+ 7 + 7 + 



1 



0-1 — 6 0-2 + 6 0-3 + 6 0-4 — 6 
{k + £)2 - Au^ ^ 



{k + iy 



{P + ke- 2cj2) (A;2 + H - 2u;2) + 2cu2) (A;2 + £2 + H - 2c<;2) 

-2(2^2)4 ^ ^^(2o;2)3 + ^^(2^2)2 ^ ^^(2a;2) + 

(£2 + H - 2cj2) (A;2 + H - 2cj2) (H + 2cu2) (fc2 + p + ki- 2u;2) ' 
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where Xj (1 ^ J ^ 4) are polynomials in k and i. For each k and £, the numerator vanishes 
at most at four values of 2ti;^. Together with the fact ( p.52| ), we have that for G (|, |), 
D does not vanish except for in a finite subset. 
The denominator in (|2.47|) has the representation: 



\U\^ _\V\^ = Re{{U + V){U -V)], 



where 



U + V 



0-4 



iba 



0"T 



52 



+ 



U ~V = D + iha 



+ higher order terms in e 
1 1 



(To 



+ 



62 



+ 



a; 



62 



_al - 62 



62 



62 



+ 



62 



+ higher order terms in e. 



Then 



\v\ 



(T4 



0-4 + 6 



Therefore, for G (|, |), the denominator in ( p.47| ), |t/|2 — |l^|2 does not vanish except for 
C(j in a finite subset. (|2.43|) - (|2.47|) give the solution to the linear system (|2.39|) - (p.42|) for 
uj G (|, I) /S, where 5* is a finite subset. 

As in 0, K{f, f) is also a bounded bilinear map: 



is>l). 



mf, fills <c\\f\\i 

We can invert the equation 

g = f + K{fJ) 
to obtain 

f = 9 + IC{9), 

where /C is of order C>(||c/||2), (s > 1). Thus, terms like J, /) and K{f,dtf - L J) 
are cubic terms in g. □ 



Remark 2.2. In this remark, we would like to make a comparison between the above nor- 
mal form transform with that in j^, and in particular to comment on why the above normal 
form transform is necessary when singular perturbation ed^f is studied. In /Q/, the linear 
operator is replaced by Lq (i.e. setting e = in L^) in constructing normal form trans- 
form. The corresponding normal form transform is given by 

K = 2u, k3{k,i) = 0, k,{k,i) = -^^, k^{k,i) = -—^^. (2.53) 
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When such a normal form transform is applied to Equation ( p. 29 ), the singular perturbation 
term ed^f will introduce the following term in the equation for g: 

ed'MfJ) (2.54) 

which is actually an unbounded bilinear operator. Therefore, we have to work with for a 
normal form transform. On the other hand, in j^, the singular perturbation ed^ is mollified 
into a bounded pseudo-differential operator ( actually a bounded Fourier multiplier) ed"^ . The 
term ( |2.54| ) is replaced by 

ed'Mf, f) 

which is of order (9(e||/||^) sufficient for the estimate on the size of the local stable manifold 
ofQ^. Although the normal form transform (|2.43| )-( p.47| ) has a more complicated expression 



than ( |2.53| ), they have the same asymptotic nature in k and i. ed\f only introduces small 
perturbations in the expressions of Kj{k,i) (1 < J < 3). 

We apply the normal form transform given by (|2.43|) - (|2.47|) to the full equation (|2.29|) , 
and the full system ( |2.27| )-( p.29 ) is transformed into: 



J 



-2a{J + uj^) + 2f3VJ + cos 9 + e7^^, (2.55) 

e = -2J-eP^L= + ni (2.56) 

gt = L,g + V,g + Af, (2.57) 
where L^, Ve, TZ^ 7^2 given in ( p.30| )-( |23^ ) with f = g + ]C{g) , and 
Af = VMg) - ^{^2 - m - iM, + K{dtf - LJ, f) + /C(/, dtf - LJ) 

= VMg) - - M2) - ^^^3 + k{vj - - ^a/'s, /) (2.58) 

+ K{f,VJ -iM2-iM^), 



where M2, M2 and A^s are given in (p:3^ ), (g) and (g3|) with / = ^ + /C(^). A/" has the 
estimate, 

\\ms-0{\J\\\9\\l + 49\\l+\\9\\% {s>l). (2.59) 
2.4 Unstable Fibers 

Under regular perturbations as in center-stable, center-unstable, and center manifolds, 
and Fenichel stable and unstable fibers persist as in the standard theory. Under the singular 
perturbation, what are the objects that persist? We start with the linear operator L^. 

2.4.1 The Spectrum of 

The spectrum of consists of only point spectrum. The eigenvalues of are: 



;i± = -e(a + P) ± k^Auj^ - P, (A; = 1,2,...). (2.60) 

When uj G (|, l), only /i^ are real, and /i^ are complex for k > 1. When uj G (l, |), only 
Hi and H2 are real, and /x^ are complex for k > 2. 
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fx 



e = 



fx 



e > 



Figure 2.3: The point spectra of the hnear operator L^. 



Remark 2.3. The main difficulty introduced by the singular perturbation ed^f is the break- 
ing of the spectral gap condition. Figure \2.^\ shows the distributions of the eigenvalues when 
e = and e 7^ 0. It clearly shows the breaking of the stable spectral gap condition. As a 
result, center and center-unstable manifolds do not necessarily persist. On the other hand, 
the unstable spectral gap condition is not broken. This gives the hope for the persistence 
of center-stable manifold. Another case of persistence can be described as follows: No- 
tice that the plane U is invariant under the PNLS flow ( |1 . 1|) . When e = 0, there 
is an unstable fibration with base points in a neighborhood of the circle S^j ( p.3| ) in U, as 
an invariant sub-fibration of the unstable Fenichel fibration with base points in the center 
manifold. When e > 0, the center manifold may not persist, but U persists, moreover, the 
unstable spectral gap condition is not broken, therefore, the unstable sub-fibration with base 
points in U may persist. This is the topics of this subsection. Since the semiflow generated 
by PNLS (|1.1|) is not a perturbation of that generated by the unperturbed NLS due to 
the singular perturbation ed^, standard results on persistence can not be applied. 

From now on, we will take the case of two unstable eigenvalues as our example to 
conduct the arguments. The case of one unstable eigenvalue is easier. The eigenfunctions 
corresponding to the real eigenvalues are: 

e± = e±^^*cosA;x, ^±^^. = ^^^^/^HE^ k = 1,2. (2.61) 

2,uj 

Notice that they are independent of e. The eigenspaces corresponding to the complex 
conjugate pairs of eigenvalues are given by: 

= span^jcos fcx}. 

and have real dimension 2. 
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2.4.2 The Set-Up of Equations 

For the goal of this subsection, we need to single out the eigen-directions (|2.61|) . Let 

±;fc=l,2 

where are real variables, and 

{h) = {hcosx) = (/icos2x) = 0. 
In terms of the coordinates (^^, J, 9, h), (|2.55|) -( p37D can be rewritten as: 

4^ = f^t^t + v.^^t + K\ {k = 1,2), 



J = e 

e = -2 J - ef3 



2a{J + cu^) + J + cos 9 
sin 9 



2 5 



7^: 



/it = L^h + K/i + A/", 
4" = -"."C + V^Cu + AT,-, (A; = 1, 2), 



(2.62) 
(2.63) 

(2.64) 

(2.65) 
(2.66) 



where //^ are given in (|2.60|) , and M are projections of M to the corresponding direc- 
tions, and 



K+a = 2c,J(4+ + C) + e/3 



sin 9 



VJ + 



k'^k 



V7C = -2c.J(4+ + C) + e/? 



sin 9 



k 



Ck 



v/4cj2 - k' 



kV4:UJ^ - k- 



k )i 



2u? 



■n Ck 



kV4:Uj'^ - k"^ 



2.4.3 Statement of the Unstable Fiber Theorem 

The main unstable fiber theorem can be stated as follows. 



Theorem 2.1. There exists an annular neighborhood A of the circle S^i (|2.3| ) in H ( p.l|) , 

for anyp G A, there is a local unstable fiber J-'^ which is a 2D surface. JF+ has the following 
properties: 

1. JF+ is a smooth surface m || ||„ norm, Vn > 1. 

2. Tp is also smooth in e, a, (3, uj, and p in || ||„ norm, for any n > 1, e E [0,eo) 
for some eo > 0. 

3. p E J-'p, Tp is tangent to span{e]^,e^} at p when e = 0, where ik = 1,2) are 
defined in ( p.61| ). 
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Figure 2.4: The bump function rj. 



4- J-'p has the exponential decay property: Let S'* be the evolution operator of (|2.62| )- 



where /i"^ = min{/i]^, yU^}. 
^- {'^p}peA forms an invariant family of unstable fibers, 
5*^;cJ-+p, VtG[-T,0], 
and WT >0 (T can be +00), such that S^p E A, Vt e [-T, 0]. 

2.4.4 Proof of the Unstable Fiber Theorem 

There are two main approaches in estabhshing invariant manifolds and fibrations: 1. 
Lyapunov- Perron's method 0, 2. Hadamard's method 0. Here we will adopt the Lyapunov- 
Perron's method, pay special attention to non-standard applications of the method, and 
focus on the difficulties generated by the singular perturbation ed^. 



Definition 1. For any S > 0, we define the annular neighborhood of the circle ( |2.^ ) as 

A{6) = {{J,9) I |J| < 6}. 

To apply the Lyapunov- Perron's method, it is standard and necessary to modify the J 
equation so that A{4:6) is overflowing invariant. Let rj G C°°{R,R) be a "bump" function: 

'0, in (-2,2) U (-CX), -6) U (6,cx)), 
1, in (3, 5), 
-1, in (-5, -3), 

as shown in Figure |2]^, \ri'\ < 2, \ri"\ < C. We modify the J equation (|2.63|) as follows: 



j = ebr]{J/6) +e -2a{J + u^) + 2(3VJ + 00^ cos 6 + e7^^, (2.67) 

where b > 2{2auj^ + 2(3uj). Then ^(45) is overflowing invariant. There are two main points 
in adopting the bump function: 
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1. One needs A{4:S) to be overflowing invariant so that a Lyapunov- Perron type integral 
equation can be set up along orbits in A{4:6) for t G {—oo, 0). 



2. One needs the vector field inside A{26) to be unchanged so that results for the mod- 
ified system can be claimed for the original system in A{S). 



Remark 2.4. Due to the singular perturbation, the real part of /i^ approaches — oo as 
k ^ oo. Thus the h equation ( 2.65|) can not be modified to give overflowing flow. This 
rules out the construction of unstable fibers with base points having general h coordinates. 



For any (70,6*0) E A{4:S), let 

J=J^{t), e = e4t), tG(-oo,0], 



(2.68) 



be the backward orbit of the modified system ( ^.671 ) and ( |2.64| ) with the initial point (Jo, ^o) 
If 

it) ,Mt) + Jit),e.it) + m , /i(t) , 4- (t) ) 

is a solution of the modified full system, then one has 



Ut 



f^t^t + K^ (fc = l,2) 
Au + F, 



where 



u 







(0 











o\ 




(Fj\ 







-2 
















Fe 


h 


, A = 
















, F = 


Fh 


er 




















F{ 


u2-y 






















F^ = v^it + K 



'k ' 



Fj = th [r/(J/5) -r/(J,(t)/5)] + e [-2aJ + 2/?VJ + cj2cos0 



- 2l3^J,{t)+uj'^ cosO^it) 



+ e7^: 



2 , 



V,h + Af, 



(2.69) 
(2.70) 



j = .Ut) + .j, e = e,{t) + ~e. 

System (|2.69|) - (p.70|) can be written in the equivalent integral equation form: 

itit) = C(io)e'^^(*-*") + r e^^(*--)F+(r)rfr, {k = 1, 2) 



(2.71) 



to 
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Jto 



(2.72) 



By virtue of the gap between /x^ and the real parts of the eigenvalues of A, one can introduce 
the following space: For a G /^^ — ^^^{f^t ^ f^t}^ ^ind n > 1, let 



a,n = U(t) = (4+(t),«(t)) 



t G (—00,0], 5'(t) is continuous 



in t in i/" norm , \\g\l^^ = supe'-^^i^ |4t(t)| + 



\u 



< oo 



t<0 



fc=l,2 



Ga,n is a Banach space under the norm || ■ ||o-,n- Let Ba^nij") denote the ball in Ga^n centered 
at the origin with radius r. Since A only has point spectrum, the spectral mapping theorem 
is valid. It is obvious that for t > 0, 

||e^*M||„ < C(l + t)||M||„, 



for some constant C. Thus, if g{t) G E^^nij), r < oo is a solution of (|2.71| )- (|2.72|) , by letting 
to — > —CO in ( ^.721 ) and setting = in ( |2.71| ), one has 



it it) = 4+(0)e^^' + f e^Ut-^F^{T)dT, {k = 1, 2) 
Jo 

u{t) = I e^^'-^^F{T)dT. 

J — oo 



(2.73) 
(2.74) 



For g{t) G B^^nij), let T{g) be the map defined by the right hand side of (|2.73| )- (|2.74|) . 
Then a solution of ( |2.73| )-( p.74| ) is a fixed point of F. For any n > 1 and e < 5^, and 
5 and r are small enough, and F are Lipschitz in g with small Lipschitz constants. 
Standard arguments of the Lyapunov-Perron's method readily imply the existence of a 
fixed point g^. of F in Bcr,n{r)- The difficulties lie in the investigation on the regularity of 
(7* with respect to {e,a,P,u, Jo,9o,C,^{0)). That is our focus. The most difficult one is the 
regularity with respect to e due to the singular perturbation, which is our further focus. 
Formally differentiating g^: in (|2.73|) - (|2.74|) with respect to e, one gets 



a.it) 



Ueit) 



=1,2 



r 



+ nt{t), (A; = 1,2) 



d^F-u,+ Y,^t^-^t 



=1,2 



(r)rfr + 7^(^), 



(2.75) 



(2.76) 



where 



Kit) = itmtX'^' + f <eit - T)e^"^^'-^F^{T)dT 

Jo 

+ f e^i^'-^\d,F^ + du,F^ ■u,^,]{T)dT, 
Jo 



(2.77) 
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+ 



(t - T)A,e^^'-^^F{T)dT 
+ 



A. 



-{a + k'^] 

/O 



0- 



\0 



1,2, 






a + dl 









(a + 1) 












-(a + 4)/ 



(J„^„ 0,0,0)^ 



(2.78) 
(2.79) 

(2.80) 

(2.81) 



where T = transpose, and (J*, 6'*) are given in ( p.68|) . The troublesome terms are the ones 
containing A^ or -u^^^ in ( |2.77| )-( pI78D . 

\\AF\\n < C\\F\\n+2 < ci\\u\\n+2 + 1^1), (2.82) 

k=l,2 

where c is small when ( ■ ) on the right hand side is small. 



dj^Fj ■ J, 



-b[v'{J/6)-v'iJ./S)]J.„ 



cos 9 



-P 



COS 6* 



+ easier terms. 



\dj,Fj-J, 



<-b sup |V'([7J* + (1-7)^1/5)1 1^1 1^*,. 

0<7<1 

+ e(3C{\J\ + \0\)\J^,e\+ easier terms 
< Ci(|J| + |^^|)|J*,e| + easier terms, 



sup e ''^^\dj,Fj ■ J^^^l < Ci sup[e 



t<0 

+ easier terms 



(|J| + |^|)]sup[e' 



uti 



where sup^^g e*^*! J=K^e| can be bounded when e is sufficiently small for any fixed u > 0, 
through a routine estimate on Equations ( p.67| ) and ( p.64| ) for {J^..{t),9^.(t))- Other terms 
involving u^^f_ can be estimated similarly. Thus, the || ||(j,n norm of terms involving u^^^ 
has to be bounded by || || cr+p,™ norms. This leads to the standard rate condition for the 
regularity of invariant manifolds. That is, the regularity is controlled by the spectral gap. 
The II ||o-,n norm of the term involving has to be bounded by || ||o-,n+2 norms. This is a 
new phenomenon caused by the singular perturbation. This problem is resolved by virtue 
of a special property of the fixed point g^, of F. Notice that if > cxi, n2 > ni, then 
Go-2,n2 C Go-i,ni- Thus by the uniqueness of the fixed point, if is the fixed point of T in 
G(T2,n2, 9* is also the fixed point of T in Gai,ni- Since g^: exists in Go-,n for an fixed n > 1 
and 0" G (iqq, ^ 10z>) where u is small enough. 



I'^fc ||(T,n < C*! + C'2||5'*||cr 
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||'^||o-,n ^ C3||5'*||cr,n+2 + C4 || fi'* || o-+!>,n + C^, 

where Cj (1 < j < 5) depend upon ||5'*(0)||„ and ||5'*(0) ||„+2- Let 

M = 2(||7^+||,,„ + ||7^|U,„), 

and r' denote the hnear map defined by the right hand sides of ( |2.75D and ( ^.761 ). Since the 
terms d^F^ , d^+F^ , d^F, and d^+F all have small || ||„ norms, T' is a contraction map on 
B{M) C L(7l, Ga,n), where B{M) is the ball of radius M. Thus T' has a unique fixed point 
Next one needs to show that g^:^^ is indeed the partial derivative of g^, with respect to 
e. That is, one needs to show 



lim 



l^^e + Ae) - g^{e) - g^^^AeW^^, 
Ae 



0. 



(2.83) 



This has to be accomplished directly from Equations (|2.73| )- (|2.74|) , (|2.75|) -( pT6D satisfied 
by g^ and g^,^^. The most troublesome estimate is still the one involving A^. First, notice 
the fact that e"^^^ is holomorphic in e when e > 0, and not differentiable at e = 0. Then, 
notice that g^ G G'o-,n for any n > 1, thus, e^^'^g^, is differentiable, up to certain order m, in 
e at e = from the right, i.e. 

{d^/dere^^"g,l=o 

exists in if". Let 

^(t, Ae) = e(^+^^)*^'f7* - e^'^'g, - {Ae)tdle''^' g. 



w{Ae) 



where t > 0, Ae > 0, and 

^Ae) = e(^^)*^^^7*-^?*-(Ae)t9^^7*. 
Since w(0) = 0, by the Mean Value Theorem, one has 

dw 



|w(Af 



where at A = 0, -4- = -j^, and 



\w{Ae)-w{mn< sup 

o<A<i "^e 

d+ 



(AAe 



lAe 



dw 



t[e^'''^'''^d%-dlgA. 



Since ;^(0) = 0, by the Mean Value Theorem again, one has 



dAe 

, dw 
'dA^ 



(AAe) 



dw 



din d'^in 

^^^') - T^^^)\\n < sup ||;^(AiAAe)|U|A||Ae| 
rfAe o<Ai<i "Ae^ 



where 



d'^w 
dAe^ 
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Therefore, one has the estimate 



71+4 • 



(2.84) 



This estimate is sufficient for handhng the estimate involving A^. The estimate involving 
ti^, e can be handled in a similar manner. For instance, let 



then 



Ae) = F{u^{t, e + Ae)) - e)) - Ae9„,F ■ m,,. 



Ae)||^.„ < |Ae|^ sup || [m,^,^ ■ c^^-F ■ + c?«.-F ■ M*,«](AAe)||<^,„ . 

0<A<1 



From the expression of F ( |2.7CI| ), one has 

11^*, e ' ^M*-^ ' ~l~ 9u:tF ■ Uif^e(\\u,n 

< Ci||5f*||^+2£;,n[(supe''*|M=^,e|)^ + sup e^""* 



t<0 



t<0 



and the term [ ] on the right hand side can be easily shown to be bounded. In conclusion, 
let 

h = g^{e + Ae) - g^{e) - g^^^Ae, 
one has the estimate 

||^||cr,n ^ '^ll^llcr,™ + |AepC'(||(y'^,||o-,n+4; || 5'* || o-+2!>,n) , 

where k is small, thus 

||^||cr,n ^ 2| AepC'(||(7=K||o-,n+4; || S'* || cr+2i>,n) • 

This implies that 

lim \ = 0, 



which is (| 

Let , 

one has 



Let 5'*(t) = (^fc^(i!^), w(t)). First, let me comment on 4" 



5+(0)=0,e=0 



0. From (IZTl) 



du 



dg* 



'9^(0)'"-'^" ■^^"94^(0)"-"' 

by letting ^+(0) ^ and e 0+, ki 0. Thus 

du 



5^(0) 



0. 



5+(0)=0,e=0 



I shall also comment on "exponential decay" property. Since 11^7*11^+ ^ 



15'* 



< re—\ Vt < 0. 
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Definition 2. Let g^{t) = {^^ (t) , u{t)) , where 

J — oo 

depends upon ^^(0). Thus 
4^(0)^ MO), 

defines a 2D surface, which we call an unstable fiber denoted by , where p = {Jo,6o) is 
the base point, ^^(0) G [— r, r] x [— r, r]. 



Let 5* denote the evolution operator of (|2.69|) - (p.7CI|) , then 
Vt < 0. 



S^J-'p C J-'stp, 



That is, {JF+}pg^(45-) is an invariant family of unstable fibers. The proof of the Unstable 



Fiber Theorem is finished. 



□ 



Remark 2.5. // one replaces the base orbit ( J*(t), by a general orbit for which only 

II ||„ norm is bounded, then the estimate (p^.82D will not be possible. The \\ ||o-n+2 norm of 



the fixed point will not be bounded either. In such case, g^, may not be smooth in e due 
to the singular perturbation. 

Remark 2.6. Smoothness of g^: in e at e = is the key point of the entire argument in 
this article. In the global theory in later sections, information is known at e = 0. This key 
point will link "e = " information to "e 7^ " studies. Only continuity in e at e = is not 
enough for the study. The beauty of the entire theory is reflected by the fact that although 
e^^^ is not holomorphic at e = 0, e'^^^g^: can be smooth at e = from the right, up to certain 
order depending upon the regularity of g^, . This is the beauty of the singular perturbation. 



2.5 Center- Stable Manifold 

We start with Equations (p:62D-([2:6^), let 



h 



h 



and let En{r) be the tubular neighborhood of 

E^{r) = {{J,9,h,^^,^,) e I \\v\\n < r}. 



(2.85) 



3D: 



(2.86) 



En{r) is of codimension 2 in the entire phase space coordinatized by {C^^, J, 9, h, C,-^ , ,^2 
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2.5.1 Statement of the Center-Stable Manifold Theorem 

Theorem 2.2. There exists a smooth codimension 2 locally invariant center-stable 
manifold W^^ in H"' for any n > 1. can he represented as the graph of a function 
C,^ : En{r) — > R'^, for some r > 0. 

1. At points in the subset VF^^4 of , is smooth in e for e G [0, eo) and some 
eo > 0. That is, if v & En+4{r) C En{r), then is smooth in t, in norm. 
Moreover, d^^^iy) is uniformly hounded in v E En+i{r) and e G [0,eo). 

2. W^^ is smooth in {a,(3,Lu). 

3. The annular neighborhood A in Theorem is included in , i. e. ^^{J, 0, 0, 0, 0) = 
0. Along the circle ( p.3|) , W^'^ is tangent to E^ir) when e = 0, i.e. i9„C^(0, 0, 0, 0) = 
when e = 0. W^'' is close to En{r), i.e. \\dy^'^\\ < Cr. 



Remark 2.7. regularity in e is crucial in locating a homoclinic orbit. As can be seen 
later, one has detailed information on certain unperturbed (i.e. e = 0) homoclinic orbit, 
which will be used in tracking candidates for a perturbed homoclinic orbit. In particular, 
Melnikov measurement will be needed. Melnikov measurement measures zeros ofO{e) signed 
distances, thus, the perturbed orbit needs to be 0{e) close to the unperturbed orbit in order 
to perform Melnikov measurement. 



2.5.2 Proof of the Center-Stable Manifold Theorem 

Let X e C°°(i?,i?) be a "cut-ofF' function: 



X 



0, in (-00, -4) U (4, oo) 

1, in (-2,2). 



We apply the cut-off 

Xs = x{\\v\\n/S)x{^t/S)x{a/S) 



to Equations (|2.62| )- (|2.66|) , so that the equations in a tubular neighborhood of the circle 
5*^^ ( |2.3| ) are unchanged, and linear outside a bigger tubular neighborhood. The modified 
equations take the form: 



vt = Av + F, 
where A is given in ( |2.7U| ), 



K = Xs[V,^^t+Afj^], 



(2.87) 
(2.88) 



F = {Fj, Fe, Fh, F^,F^f, T = transpose. 



Fj = X5 e 



-2a{J + o;2) + 2/3VJ + a;2 cos ^ + 7^ 
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Fe = Xs 



^ sin 9 



VJ + 



FH = Xs[V.h + U], (A; = 1,2) 
Fu=X6[Vak+Kl 
Equations (g^STD-f^) can be written in the equivalent integral equation form: 



to 



v{t) = e 



A{t~to) 



(to)+ fe^'^'~^^F{r)dT. 

Jto 



(2.89) 
(2.90) 



We introduce the following space: For a G /^^ = mm{fif, /x^}, and n > 1, let 



G'.,.= ^(t) = (4+(t),i;(t)) 



t G [0, oo), g{t) is continuous in t 



in i7"norm, \\g\\,,n = supe~'^*[^ 1^(^)1 



V{t n < OO 



t>0 



k=l,2 



Ga,n is a Banach space under the norm || ■ ||o-,n- Let Aa,n{r) denote the closed tubular 
neighborhood of S^^ (|2.3|): 



snpe-'lY^ l4t(t)| + ||^(t)IU]<r 



k=l,2 



where w is defined in (|2.85|) . If g{t) G Aa,n{r), r < cx), is a solution of (|2.89|) - (|2.90|) , by 
letting to +00 in ( ^.89| ) and setting to = in ( ^.90| ), one has 



&)= f e'^i^(*--)F+(r)rfr, (A; = 1,2) 

J +00 

v(t) = e^\;(0) + [\^^'~^^F{T)dT. 



(2.91) 
(2.92) 



For any g(t) G ylo-,n('^), let T{g) be the map defined by the right hand side of ( |2.91| )- (|2^9^) . 
In contrast to the map F defined in (|2.73|) - (|2.74|) , F contains constant terms of order 0{e), 
e.g. Fj and Fg both contain such terms. Also, Aa,n{i^) is a tubular neighborhood of the 
circle 5"^^ ( p.3|) instead of the ball B^^ni^) for F. Fortunately, these facts will not create any 
difficulty in showing F is a contraction on Acr,n{f)- For any n > 1 and e < 5^, and 6 and 
r are small enough, F^ and F are Lipschitz in g with small Lipschitz constants. F has a 
unique fixed point g^, in Aa,n{r), following from standard arguments. For the regularity of 
g^ with respect to (e, a, P, co, f (0)), the most difficult one is of course with respect to e due 
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to the singular perturbation. Formally differentiating in (p.91|) -( p^92D with respect to e, 
one gets 



vit) 



+00 



=1,2 



(r)c/r 



+ Rtit), ik = l,2) 



=1,2 



(r)cir + i?(t), 



(2.93) 



(2.94) 



where 



m 



f f^tit-T)e^^^'-^F,^iT)dT+ f e'^^(*--)9,F+(r)rfr, 
J +00 J +00 

M,e%(0) + fit - r)A,e^(*~")F(r)rfr + f e^(*-")9,F(r)rfr, 
Jo Jo 



(2.95) 
(2.96) 



and n'l^ and are given in ( |2.79| )-( p.81| ). The troublesome terms are the ones containing 
Af in i^M)- These terms can be handled in the same way as in the Proof of the Unstable 
Fiber Theorem. The crucial fact utilized is that if f (0) G H"'^, then is the unique fixed 
point of r in both Go-,ni and Ga^n2 for any n2 < rii. 

Remark 2.8. In the Proof of the Unstable Fiber Theorem, the arbitrary initial data in 
( p.73| )-( p.74 ) are ^^(0) (A; = 1,2) which are scalars. Here the arbitrary initial datum in 
( [^.91| )-( p.92D is f (0) which is a function of x. //f (0) G i/"^ but not iJ"^ for some ni > n2, 
then g^ ^ Ga,ni, in contrast to the case of (|2.73|) - (|2l7^) where g^ G G„^n for any fixed 
n>l. The center-stable manifold W^^ stated in the Center- Stable Manifold Theorem will 
be defined through f(0). This already illustrates why has the regularity in e as stated 
in the theorem. 

We have 

||-Rfe"||(T,n ^ C*!) 
||-R||o-,n < C'2||^*||(T,n+2 + C*3) 

for G Aan+2ir), where Cj (j = 1,2,3) are constants depending in particular upon the 
cut-off in F^ and F . Let F' denote the linear map defined by the right hand sides of 



( p.93| )- (|2.94| ). If f (0) G -f/'"+^ and g^ G X,n+2(r), standard argument shows that F' is 
a contraction map on a closed ball in LiR,Ga,n)- Thus F' has a unique fixed point g^^e- 
Furthermore, if f(0) G and g^.. G Aa,n+4ir), one has that ^^.^^ is indeed the derivative 

of ^* in e, following the same argument as in the Proof of the Unstable Fiber Theorem. 
Here one may be able to replace the requirement f (0) G and g^: G Aa,n+4ir) by just 

f (0) G and g^: G Aa,n+2ir)- But we are not interested in sharper results, and the 

current result is sufficient for our purpose. 
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Definition 3. For any v{0) G En{r) where r is sufficiently small and En{r) is defined in 
( p.86| ), let g^:(t) = {C,^{t),v{t)) be the fixed point ofV in Ga^n, where one has 



f e^^^'-^^F+{T)dT, (A; = 1,2) 



which depend upon v{0). Thus 

^(0)^4^(0), (A: = 1,2) 
defines a codimension 2 surface, which we call center-stable manifold denoted by Wf^ . 

The regularity of the fixed point g^, immediately implies the regularity of Wf^ . We 
have sketched the proof of the most difficult regularity, i.e. with respect to e. Uniform 
boundedness of 9^^^ in t>(0) G i?„+4(r) and e G [0, eg), is obvious. Other parts of the detailed 
proof is completely standard. We have that W!^ is a locally invariant submanifold which 
is in (a, /5, lS). Wf^ is in e at point in the subset W^^^. From Equation ( p.91|) . Claim 3 



in the Theorem immediately follows. □ 

Remark 2.9. Let 5** denote the evolution operator of the perturbed nonlinear Schrddinger 
equation (|1 . 1| ) . The proofs of the Unstable Fiber Theorem and the Center-Stable Manifold 
Theorem also imply the following: is a map on if" for any fixed t > n > 1. 5** 
is also in {a,(3,uj). S** is in e as a map from H""^^ to H"' for any fixed n > 1, 
e G [0,eo), eo > 0. 

2.6 Stable Manifold of Q, 

As mentioned earlier, the homoclinic orbit to be located will be asymptotic to the saddle 
( |2.10|) . Dynamics on the invariant plane 11 ( p.l|) on which lives, is governed by 
Equations (|2]|)-(|2]|) which are equivalent to Equation (|2]23)-(g]5|) with / = 0. The 
eigenvalues of are given by (|2.11|) on 11 and ( p.6(]| ) off 11. Thus has three unstable 



eigenvalues of two scales: One unstable eigenvalue of order 0{^/e) with eigen-direction in 
n, the other two unstable eigenvalues of order 0{1) with eigen-directions off 11. On 11, 
has the unstable curve 0^ with approximate representation ( [2.221 ). Thus the 3D unstable 
manifold of Q^, W^{Q^) has the representation 

where J^^ is the unstable fiber given in Theorem |2.1|. The scales of the stable eigenvalues 



of Qe range from 0{t) to O{oo). The stable eigenvalue with eigen-direction in 11 has order 
(9(y^). On n, has the stable curve 0^ with approximate representation ( [^.22| ). From 
the standard stable manifold theorem, has a stable manifold Wn{Qe) in if" for any 
n > 1. In fact, the codimension 3 stable manifold of Q^, W^{Qe) intersects 11 along 0^. In 
order to locate a homoclinic orbit, we need the size of W^{Q^) large enough. Along 0^, 
the size of W^{Qe) is 0{1) sufficient for our purpose. One can view W^{Qf) as a wall with 
base 0y^. As can be seen later in the Second Measurement, one needs the size of W^iQ^) off 
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n to be of order 0{e'^), /t < 1 in order to overcome the order 0{e) "fuzz" between certain 
perturbed and unperturbed (e = 0) orbits to locate a perturbed homoclinic orbit. Starting 
from the system ( |2.27| )-( pl29D , one can only get the size of W^{Qf:) off 11 to be 0(e) from 
standard stable manifold theorems. As discussed previously in the subsection on Normal 
Form Transformation, an estimate of order 0{e'^), k < 1 can be achieved if the quadratic 
term N'2 (|2.34 ) in ( p. 29 ) can be removed through a normal form transformation. Such a 



normal form transformation has be found in that subsection. 

Theorem 2.3. The size of W^{Q,) off U is of order 0{^) for G (|, |) /S, wh ere S is 
a finite subset. 

Proof. For uj G |) /S, where is a finite subset, we apply the normal form transform 
given by (|]43])-(|]4^) to Equation ( ^1291) , then the system (|]23)-(|2]2|) is transformed 



into the system ( p.55| )-( p.57| ). By virtue of the estimate ( |2.59| ), the theorem follows from 



standard argument. For details, see y]. □ 

As discussed in the subsection on Center-Stable Manifold, the center-stable manifold 
Wf^'^ is unique. Thus W,^{Q^) is a codimension 1 submanifold of W^^. 
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3 Global Theory 



Global Theory is referred to a theory global in phase space, which includes integrable 
theory, Melnikov measurement, and the so called second measurements. These are tools 
necessary in locating a homoclinic orbit. 

The entire process of locating the homoclinic orbit can be briefly summarized as follows: 
The integrable theory will provide explicit representations for certain family of homoclinic 
orbits asymptotic to periodic orbits on the invariant plane 11. Local unstable fiber theorem 
will provide ways of picking orbits in the local unstable manifold of Q^-, that are close 
to certain unperturbed homoclinic orbits. Our main strategy is to use the unperturbed 
homoclinic orbits to trace the candidates for a perturbed homoclinic orbit. The procedure 
is splinted into two steps: 

Step 1. Find an orbit that is in W{Q,) n W^' . 

Step 2. Find out when this orbit is also in W^iQ^:) nH^^((5e), where W^{Q^) is a codimension 
1 submanifold of W^^ . 

Step 1 will be accomplished through Melnikov measurement. The Melnikov vectors will 
be provided by integrable theory. The Melnikov integrals will be evaluated along the 
unperturbed homoclinic orbits mentioned above. In contrast to the work ||^, the new 
feature in Step 1 is that W^^ is not in e everywhere rather only at its subset W^^^. This 
difficulty is overcome by the fact that C for any fixed n > 1 by virtue of the 

unstable fiber theorem. Step 2 will be accomplished by the so called second measurement. 
It turns out that one can trace the perturbed orbit in IV" (Qe)niy^'* through an unperturbed 
homoclinic orbit to an order 0{e\ lne|) neighborhood of 11 In order to check when 

this orbit can be in W'^{Q^) fl W^{Q^), one needs the size of W^{Q^) off 11 to be large 
enough, and 0{^/e) is sufficient. 

3.1 Integrable Theory 

Consider the integrable ID cubic focusing nonlinear Schrodinger equation (e = in ( |1 . 1| ) ) , 



iqt = Qxx + 2[|gp - uj'^]q . 



(3.1) 



Its Lax pair is given by the Zakharov-Shabat linear system. 



(3.2) 
(3.3) 



where 
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3.1.1 Isospectral Theory 

Focusing one's attention on the spatial part (|3.2| ) of the Lax pair ( p.2| , |3.3|) , one can define 
the fundamental matrix solution M{x), s.t. M(0) is the 2x2 identity matrix. Then the 
Floquet discriminant A is defined as 

A = trace M{2n) . 

A = A(A, g), as a functional in q for any A G C, provides enough functionally independent 
constants of motion to make NLS integrable in the classical Liouville sense. For each 
fixed q, there is a sequence of special points {A^, j E Z} of X E C called simple points for 
which |A(A^,g)| = 2. There is also a sequence of critical points {A^, j G of A G C for 
which |^A(A^^,g) = 0. When some A^ coincides with some A^, a double point is formed. 
The geometric multiplicity is the dimension of the eigenspace of ( p. 21) at the double point. 

Definition 4. The sequence of constants of motion Fj is defined as 

F, = A(A,^,g) , jeZ . (3.4) 

Fj's provide a sequence of Melnikov functions. More importantly, the gradients of F/s, 
which will be the Melnikov vectors, have a simple representation, 

'^^-^ f ^^^^ ^ atA^A?, jeZ, (3.5) 



6q w{ij+,ij-) V -i^Ii^'i 

where q = {q,q)^, tp"^ = {ipf,tp2)'^ ^"^0 eigenfunctions at A = A^, and W{^p~^ ,ip~) = 
ipfip2 ~ "^ti^i is Wronskian. For more details on the isospectral theory of NLS, we 
refer the readers to [^. 

3.1.2 Backlund-Darboux Transformation 

The particular form of the Backlund-Darboux transformation for NLS ( p.l|) , that is useful 
for our purpose, is due to David Sattinger and V. Zurkowski 0. 

Theorem 3.1. Let q{t,x) be a solution of NLS u is a complex double point of ge- 

ometric multiplicity 2. Let 0^ be two linearly independent eigenfunctions of the Lax pair 



( \3. 3i , \3. 5| j at \ = V. Denote by cj) the general solution 

= 0(t, X, Z/, C+, C_) = C+0+ + C_0~ , 

We use ip to define a Gauge transformation matrix 

G = G(A;i.,0)=r(^ ^"'^ ^^-^^ 

where 

^ \ -02 
^2 01 
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Then we define Q and \1/ by 



Q = q + 2{v-v) 



>W2 



l0lP+l02P 



and 



where ip solves the Lax pair (^~^\3. 3^ ) at (X,q)- Then \1/ solves the Lax pair ^3. 2{ , \3. Sj ) at 



(X, Q ), and Q also solves NLS ( \3. j|j 



3.1.3 Figure Eight Structures 

Consider the special solution of NLS ( |3.1|) , 

= ae'^W , e{t) = -[2{a^ - uj^)t + 7] . 
The corresponding Floquet discriminant is given by 

A(A,gc) = 2cos(27rA;) , = Va^ + , 
and two eigenfunctions (Bloch functions) are 



(3.7) 



ae 2 



(±A;- A)e-*^ 



exp{±i2Xkt ± ikx} . 



(3i 



When k is real, to have temporal growth (and decay) in ip^, one needs A to be purely imag- 
inary. The temporal growth (and decay) in is connected to the linear instability of qc, 
since quadratic products of solve linearized NLS [Q. The temporal growth is also neces- 
sary for constructing homoclinic solutions through the Backlund-Darboux transformation. 
Specifically, the double points of A are given by 



k = Va2 + A2 = j/2 , J G Z/{0} . 
If one requires that a lies in the interval 
a G (1/2,1) , 

then there is only one pair of complex double points 
X = ±u = ±ia , a = a/o^ - 1/4 . 

If one requires that a lies in the interval 
a G (1,3/2) , 

then there are two pairs of complex double points 



A = ±1' = ±ia , and A = ±/> = ±ia , a = V 0? — I ■ 
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Figure 3.1: Figure eight structure of noneven data with one unstable mode. 



Next we will construct homoclinic orbits, starting from the special solution Qc, through the 
Backlund-Darboux transformation. Notice that building the Backlund-Darboux transfor- 
mation at A = v.s. at A = and at A = z> v.s. at A = — z> always lead to equivalent 
results. We will choose X = u and X = 0. 



One Pair of Complex Double Points Case 

Let 0^ = ip'^itjX, p) defined in (^.81), and let 



(3.9) 



Applying the Backlund-Darboux transformation given in Theorem |3.1| , one gets a new 
solution, 



Q = gc 1 + sin {}o sechr cos y 



cos 2i}o — i sin 2'dn tanh r 



sin {}q sechr cos y 



where 



c+/c- = e^+*'^ , ^ + ^ = ci&"^° , T = 2at- p , y = x + d - + n /2 



As t — > ±CXD, 



(3.10) 



(3.11) 



(3.12) 



Thus Q is asymptotic to qc up to phase shifts as t ^ ±oo. We say Q is a homoclinic orbit 
asymptotic to the periodic orbit given by q^. For a fixed amplitude a of qc, the phase 7 of 
qc and the Backlund parameters p and ?9 parametrize a 3-dimensional submanifold with a 
figure eight structure. For an illustration, see Figure If one restricts the Backlund 

parameter ^9 by "i? — "i^o + 7r/2 = 0, or vr, one gets Q to be even in x, 



Q 



1 ± sin {}o sechr cos x 



cos 2{}q — i sin 2'do tanh r =1= sin {Iq sechr cos x 



(3.13) 
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Figure 3.2: Figure eight structure of even data with one unstable mode. 



where the upper sign corresponds to 0. Then for a fixed amphtude a of Qc, the phase 7 of 
Qc and the Backlund parameter p parametrize a 2-dimensional submanifold with a figure 



eight structure. For an illustration, see Figure 3.2 



Two Pairs of Complex Double Points Case 



Let 0^ = ip'^itjX, z>) defined in (p.8|) , and let 



C+0+ + C' 



(3.14) 



In this "two pairs of complex double points" case, to get the complete fohation of the 
figure eight structure, one needs to iterate the Backlund-Darboux transformation. First 
one needs to apply the Backlund-Darboux transformation at X = u, then one needs to 
iterate the Backlund-Darboux transformation at A = z/. Switching the order between u 
and z> leads to the same result. At A = z/, the Gauge transform G = G{X; u, 0) ( |3.6| ), then 
one defines 



± 



Let 



•I = G{u; u, 



(3.15) 



(3.16) 



After an iteration on the Backlund-Darboux transformation, one gets the solution of NLS 
( ^^1) with the representation. 



Q = q^ + 2{u- u) 



+ 2(z>-z>)- 



$1$2 



Explicit formula for Q is, 
W2 sin da 



Q = Q + qc 



where Q is given in ( p.lOp 



Wi 



(3.17) 



(sin'i9o)^(l + sin-i^o sechrcosy)^ H — (sin 2^9o)^(sechr)^(l — cos2y) 
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[1 + sin^^o sechfcosy) 



— - sin 2{}q sin 2^0 sechr sech'f(l + sin -(^o sechr cos y) sin y sin y 



+ (sin^^o 



1 + 2sin'(9o sechr cosy + [(cosy)^ — (cos -i^o)^] (sechr)' 



[1 + sin-i^o sechr cosy) 



— 2 sin -i^o sin i^o 



cos cos t^o tanh f tanh r + (sin {}() + sechr cos y) 



(sin do + sechr cos y) 



[1 + sin-i^o sechr cosy) 



W2 = 



+ 



— 2(sin'(9o)^(l + sin^^o sechr cosy) ^ + -(sin 2790)^ (sechr)^(l — cos2y) 



• (sin do + sechr cos y + i cos do tanh f) 
+ 2(sin'i9o)^(— cos-i^o tanhr + i sini^o + ^ sechr cos y)^ 
■(sin'(9o+ sechr cos y — z cos -i^o tanh f) 
2 sin '(9o(sin i)o + sechr cos y + i cos i^o tanh r) 

2 sin ( 1 + sin -(^o sechr cos y ) ( 1 + 



sin-i^o sechr cosy) 



— sin 2i}o cos 'do sechr sechr sin y sin y 



and the notations are given by 



1 + i> = ae'^° , c+/c- = 6'^+^'' , f = 4(jt - p , y = 2x + d - do + 7r/2 . 
The asymptotic phase of Q is as follows, as t ±00, 



Ti2(,3o+i?o) 



(3.18) 



Thus Q is asymptotic to qc up to phase shifts as t — > ±00. We say Q is a homoclinic orbit 
asymptotic to the periodic orbit given by qc. For a fixed amplitude a of qc, the phase 7 of 
qc and the Backlund parameters p, p, and parametrize a 5-dimensional submanifold 
with a figure eight structure. For an illustration, see Figure p.3| . If one put restrictions on 
the Backlund parameters and d, s.t. 



^ -^0 + 71/2 = < 







TT 



^9 - + 7r/2 = , 
^9-^o + 7r/2 = 7r , 
^9 - 1^0 + vr/2 = , 
d - do + 7^/2 = IT , 



(3.19) 



then Q is even in x. Then for a fixed amplitude a of qc, the phase 7 of qc and the Backlund 
parameters p and p parametrize a 3-dimensional submanifold with a figure eight structure. 



For an illustration, see Figure 3.4 
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Figure 3.3: Figure eight structure of noneven data with two unstable modes. 



Figure 3.4: Figure eight structure of even data with two unstable modes. 



3.1.4 Melnikov Vectors 

Notice that ( p.5|) evaluated at {u, Q) and (z/, Q) are linearly dependent. Same is true for 
{v.Q) or (i>,Q). 

One Pair of Complex Double Points Case 

In this case, the Melnikov vector is (|3.5|) at A = z^, evaluated along the homoclinic 
orbit Q ( |3Jq ) or (|37[^ ). 

_ . v/A2(z/)-4 / $+$2 \ 
where (cf: (PT^)), 

$± = G(z/;i/,0)0± (3.20) 

^-W,r)^(_|). 

By L'Hospital's rule, 

VA2 -4 _ v^A(z/)A"(z/) 

~ (z/-z/)Vr(0+,0-) ■ 

</) (|3.9| ) can be rewritten as 

0^ = 2Vc+c-ae^^/\i , 02 = 2Vc+c-ae"^^/2M2 , (3.21) 
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where 



Ui = cosh — cos z ~ i sinh — sin z , 
2 2 ' 

T T 

U2 = — sinh — cos(2; — i^o) + ^ cosh — sin(2; — -^o) 



where 



z = x/2 + ^/2 



and other notations have been defined in ( p.7| , p.l7|) . Finally, one gets the explicit 
representation for the Melnikov vector, 



Sq 4 



(3.22) 



Two Pairs of Complex Double Points Case 

In this case, the Melnikov vectors are ^ and ^ ( p. 5] ) at A = z/ and A = z/ respectively, 

evaluated along the homoclinic orbit Q ( |3.17D or ( |3.19| ). We know that $ is defined in 
( p. IE ). Then we use $ to define a Gauge matrix G(A; u, $). Let 

= G(z/;z>,<l)$^ , <|)(^'*) = G'(z>;z>,<l>)<l>^ , (3.23) 

(3.24) 
(3.25) 



where $^ and $^ are defined in ( p.20| , |3.15| ). Then the Melnikov vectors are 
SFi _ . v^A2(zy) -4 



p^(<|){+.*)^ $(-,*)) 



SFo 



/A2(i)) -4 



(5(f iy(${ + .*), $(-:*)) 

By L'Hospital's rule. 



iy($(+,*), $(-■*)) (zy - zy)(zy - z})(zy - u)W {(p^ , 0") 



VA2(^>)-4 



v/A(^)A^ 



iy($(+,*), $(-,*)) (z) - z>)(z} - zy)(z} - Z7)iy ((/)+, 0-) 

We know that (|3.9| ) can be rewritten as ( p.21|) . (|3.14|) can also be rewritten as 

-ie/2 



01 = 2V c+c~ae*^/^f 1 , 02 = 2V c+c~ae "'''V2 , 



(3.26) 



where 



T T 

Vi = cosh — COS z — i sinh — sin i 
2 2 
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V2 = — sinh — cos(z — -t^o) + ^ cosh — sin(;z — -(^o) 



where 

z = x + ^/2 , 

and other notations have been defined in (p.7| , |3.11| , |3.17| ). Using (|3.21| , |3.26|) , one can get 
the representation for $ ( |3.1(j| ), 

where Vi and V2 are defined as 

Vi = 1^ 'I 1^ [(z> - z/)|uip + (P - z/)|u2nt;i + (z/ - z/)MiIZ^t;2 

^2 = 1 TT— n {u - u)WU2Vl + [{l) - h')\ui\^ + {0 ~ U)\U2\^]V2 



Finally, one gets the explicit representations 



5Fi 



6F2 
Sq 



q-cSl 
-qcSf 



2 



-a-H{u - l){v -v){y- z7)v/A(z>)A"(z>) 



2 / ' 



(3.27) 
(3.28) 

(3.29) 
(3.30) 



where Si and Si {I = 1,2) are independent of the phase 7 of qc, and have the representations 

1 



^1 



^2 



^1 
S2 



{\u^\' + \u2mv,\^ + 

-(z> - u)ViV2Ul 

1 



[(i, - 0)\V,\' + (u - U)\V2\']U^ 



(3.31) 



(^-Z>)Vl\/2M^-[(z/-^)|\/lP 



(3.32) 
(3.33) 
(3.34) 



3.2 Melnikov Analysis 

Let p be any point on 0^ (|2.22|) which is the unstable curve of in 11 ( p.l| ). Let ge(0) 
and go(0) be any two points on the unstable fibers JF+ |^ and JF+ |^=o; with the same 
coordinates. By the Unstable Fiber Theorem, JF+ is in e for e G [0, cq), Cq > 0, thus 

l|ge(0)-go(0)||„+8<Ce. 
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The key point here is that JF+ c for any fixed s > 1. By Remark 2^, the evolution 
operator of the perturbed NLS equation ( |1.1|) 5* is in e as a map from to for 

any fixed n > 1, e G [0, eo), eo > 0. Also is a map on iJ" for any fixed t > 0, n > 1. 
Thus 

||g,(r) - go(T)|U+4 = WS'^iqM) - ^^(go(0))|U+4 < C,e, 
where T > is large enough so that 

go(r) e w:i, |,=o . 

Our goal is to determine when qe{T) G W^^ through Melnikov measurement. Let qe{T) and 
qo(T) have the coordinate expressions 

q,iT) = iC'%Ve), go(T) = (C'°,t;o). (3.35) 
Let (jeiT) be the unique point on which has the same t>-coordinate as qeiT), 

By the Center-Stable Manifold Theorem, at points in the subset W^^^4, W^^ is smooth 
in e for e G [0, eo), eo > 0, thus 

MT) - q,{T)\\n < C2e. (3.36) 

Also our goal now is to determine when the signed distances 

(fc = l,2) 

are zero through Melnikov measurement. Equivalently, one can define the signed distances 

4 = (VFfc(go(T)),g,(T)-g,(T)) 

^ d,Fk{qo{T)){q,{T) - UT)) 

+ dgFMT)){q,{T) - UT))-, k = l,2, 

where Fk and VFk are given in the subsection on Integrable Theory, qo{t) is the homoclinic 
orbit also given in the same subsection. In fact, q^it), q^it), qo(t) G H"^, for any fixed n > 1. 
The rest of the derivation for Melnikov integrals is completely standard. For details, see 

§ §■ 

4 = eMfc + o(e), A; = 1,2, (3.37) 

where 

kcx) /•2tt 

^2, 



Mk= / [d.Fkiqomd'Mt) - c^qoit) + (3) 

J~oo Jo 

+ dgFkiqoit)){dl'^ - + (3)]dxdt, 

where ^0(^)5 dgF^, and dqFk are given in the subsection on Integrable Theory. 
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Theorem 3.2. There exists eo > 0, such that for any e G (0,eo), there exists a domain 
T)^ C IR+ X ]R+ X R+ where uo G f)/*?, S is a finite subset, and auj < (3. For any 
{a,(3,uj) G V^, there exists another orbit in W'^{Q^) fl PV^* other than the unstable curve 
0^ ^.22^) of Q^, for the perturbed nonlinear Schrddinger equation ( |1 . 1|) . 



Proof. This theorem follows immediately from the explicit computation in the subsection on 
Evaluation of Melnikov Integrals and Second Distance, and the implicit function theorem. 

□ 



3.3 The Second Measurement 

The second measurement starts with the orbit obtained in Theorem |3.2| , i.e. q^{t) where 
qe(T) = qeiT). The goal is to determine when q^{t) is also in W^{Q^). Recall that W^{Q^) 
can be visualized as a co dimension-one wall in M/^* with base curve in 11 and with 0{^/e) 
height. Thus we have to continue to follow qe{t) and qQ{t) to a smaller neighborhood of 11. 
From the explicit expression of qo{t), we know that q'o(^) approaches 11 at the rate (9(e~^*), 



/i = min{V4u;2 _ 4Vcj2 _ 1} (3.38) 
(cf: <^M))- Thus 

distance |go(T + -I lne|),nl <Ce. (3.39) 



Lemma 3.1. For all t G 



T,T + i|lne| 



\\qe{t)-qG{t)\\n<Cie\\ne\\ (3.40) 
where Ci = Ci(T). 

Proof We start with the system (^]91])-(g]92|). Let 

q.{t) = {C''{t),r{t),e^{t),h^{t),^-''{t)), 

qo{t) = {C\t),j\t),9'{t),h'{t),^~^\t)). 
Let Ti(> T) be a time such that 

\\qe{t)~qomn<C2e\\ne\\ (3.41) 

for all t G [T, Ti], where C2 = C2iT) is independent of e. From (|3.36|) , such a Ti exists. 
The proof will be completed through a continuation argument. For t G [T, Ti], 



E + \Cu\t)\) + \\h'{t)\U < C,re-'^^^'-^\ \J%t)\ < a 



4Ve, 



k=1.2 



\J'it)\ < \j\t)\ + \J'it) - J\t)\ < \j\t)\+C2e\\ne\^ < C.^fe, (3.42) 



k=l,2 
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where r is small. Since actually qe{t), qo{t) G H"^ for any fixed n > 1, by Theorem p.2| , 

iC^'it) - C'\t)\ < C4v,{t) - voimn + Cre, (3.43) 

whenever vo(t) G En+A^r), where v^iT) = and vo(T) = vo are defined in ( |3.35|) . 

Thus we only need to estimate \\v^(t) — fo(t)||n- From ( p.92|) , we have for t G [T, Ti] that 

v{t) = e^(*-^)t;(T) + e^(*-")F(r)rfr. (3.44) 

Let Av{t) = v,{t) - vo{t). Then 

Av{t) = [e^(*-^) - e^l-°(*-^)]t;o(r) + e^^'-^^Av{T) 

By the condition ( |3.42| ), we have for t G [T, Ti] that 

||F(t) -F(t)|,=o||n < [C8v^+C9re-^'^(*-^)]e|lnep. (3.46) 

Then 

||At;(t)||„ < Cioe(t - T) + Ciire| Inep + Ci2v^(t - T)2e| lne|^ (3.47) 
Thus by the continuation argument, for t G [T, T + i| lne|], there is a constant Ci = C*i(T), 

||Af(t)||„ < Ciejlnel^ (3.48) 

□ 

By Lemma |3.1| and estimate (|3.39|) , 



distance |ge(r + ^| lne|), n| < Ce\\ne\'^. (3.49) 

Recall the fish-like singular level set given by Ti (p.20|) , the width of the fish is of order 
0{^/e), and the length of the fish is of order 0{1). Notice also that qo{t) has a phase shift 

e^, = e%T + -\\ne\)-e%o). (3.50) 

For fixed /3, changing a can induce 0{1) change in the length of the fish, 0{y/e) change in 
61, and 0(1) change in 6'°(0). See Figure for an illustration. The leading order signed 
distance from qe{T + -| lne|) to W,^{Q^) can be defined as 

= 2u [aiu9^^ + P[sm 9%0) - sin(0°(O) + 0°)]] , 

where Ti is given in ( |2.2CI| ). The common zero of ( p.37| ) and d and the implicit func- 
tion theorem imply the existence of a homoclinic orbit asymptotic to Q^. Much detailed 
arguments have been given to the signed distance d and the second measurement in and 
i- 
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Figure 3.5: The second measurement. 



3.4 Evaluation of Melnikov Integrals and Second Distance 



It turns out that to the leading order, one can evaluate ( |3.37D at go(^) where a = u. 
Our goal in this subsection is to find the common zero of ( |3.37| ) and d ( |3.51| ) . 



One Pair of Complex Double Points Case 

Ml = and d = lead to 

Ml = M^^) + aM^^) + /3cos7M(^) = , 



(3.52) 



P cos 7 



au!{A'~f) 
2sin^ 



(3.53) 



where A7 - 



-Mo, M(J) = M(-')(cu), (j = 1,2,3), and 

hoc /•2tt 



/-hOO /"ZTT 
/ d^^il' + \u2\^)-^[uldlP - uldlP]dxdt 
■00 Jo 



00 l'2-K 

2 / / /!„, |2 



00 -'0 



(l^iil + W2\ )~ \uiP - u{P]dxdt 



-00 ;'27r 

/ / (ImiI^ + |-U2p)~^['Ui — 'U2]c?X(it 

-00 Jq 



and P is given by Q = gc-P, and Q is given in ( |3.13|) . Equations (|3.52| ) and ( |3.53|) define a 
CO dimension-one surface in the space of (a,/5, tu), given by 



a 



k{uj) 



where 



k{uj) = -[2M(2) sin^ + M(3)cj(A7)][2M« ^^"^T"^"' 
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Figure 3.6: The curve of k^u). 



and its graph is plotted in Figure 



Two Pairs of Complex Double Points Case 

Mj = (j = 1, 2) and d = lead to 



M, = M. 



(1) 



+ aMf + 13 cos ^M]-"' + {3 sin 7M] 



r(3) 



(4) 



0, (j = l,2) 



(3.54) 



/? cos 7 



(4) 



2sin^ 



(3.55) 



where A7 = -4(?9o + ^o), = (cj,Ap), (j = 1,2, / = 1,2,3,4), Ap = 2aa-^p-p, 
and 



lUO 



+00 /•27r 



[^a'a^^P - SidtP]dxdt 



00 JO 

+00 l'2lT 



[SlP - SlP]dxdt 



00 JO 
+00 /•27r 



[^2^ - Sl]dxdt 



-00 JO 

00 /■27r 



00 Jo 



[5*2 + Sl]dxdt , 



and P is given hy Q = QcP, and Q is given in (|3T^ ) and (|3l9| ). iVff can be obtained from 
m[''^ {I = 1,2,3,4) by replacing Sm by Sm {m = 1,2). Equations (p.54|) and ( p. 551) define 
a codimension-one surface in the space of {a,P,uj), given by 



a 



X(^, Ap) 
P = /3(^,Ap) 



(acjA7)2(2sin^)"2 



+ (Mf V'fMf ^ + a(M}"^ + Mi^"^cuA7(sin 



(2) , A AS), 
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Figure 3.7: The surface of x(c<j, Ap). 



where 



+a;A^(2 sin ^)-i(Mf)Mf) - Mf Mf^; 



and its graph is plotted in Figure 3.7. 



3.5 Statement of the main Theorem 

Theorem 3.3 (Main Theorem). There exists a eo > 0, such that for any e e (0,eo), 
there exists a codimension 1 surface in the space of a;) G M"*" x R"*" x M"*" where 
uj G (|, '^)/ S, S is a finite subset, and au < 13. For any u;) on the codimension- one 
surface, the perturbed nonlinear Schrddinger equation ( |1 . 1|) possesses a homoclinic orbit 
asymptotic to the saddle ( |2.1CI| ). The codimension 1 surface has the approximate repre- 
sentation given in the subsection on Evaluation of Melnikov Integrals and Second Measure- 
ment. 

Proof. From the exphcit computation in last subsection and the implicit function theorem, 
dk and d are zero for the parameter values specified in the theorem . □ 
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